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Abstract 

In this paper we propose a solution for the time evolution of the island density with irreversible 
aggregation and a time dependent input of particle in the space dimensions d = 1,2. For this 
purpose we use the rate equation resulting from a generalized mean field approach. A well-known 
technique for growing surfaces at the atomic scale is molecular beam epitaxy (MBE). Another 
approach is the pulsed laser deposition method (PLD). The main difference between MBE and 
PLD is that in the case of MBE we have a continuous rate of deposition F of adatoms on the 
surface whereas in the case of PLD the adatoms are deposited during a pulse of a laser which is 
very short in comparison to the time span T between the pulses. The generalized mean field theory 
is a useful model for both MBE and PLD with the most simple approximation, point-like island. 
We show that the parameter T distinguishes the MBE regime from the PLD regime. We solve 
the rate equation for the PLD regime. We consider the time evolution of the density of immobile 
islands. For large time f > T, the PLD regime dominates the MBE regime and we find that the 
density of immobile islands grows as i 1 / 2 whereas for MBE we find the known behavior of the 
density, t 1 / 3 for d = 2 and i 1 / 4 for d = 1. We illustrate this result with Monte-Carlo simulations 
for d = 1,2. The author recognizes that in real experiments , some deviations from this simple 
point island approximation in the rate equations could arise. 
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I. INTRODUCTION 



Surface growth by pulsed deposition plays an important role in the fabrication of thin 
films. In this technique the material is ablated by a pulsed laser and then deposited in pulses 
so that bunches of many atoms arrive at the surface simultaneously^. 

Pulsed laser deposition plays an important role in various applications including the 
growth of ultra-hard carbon films, artificially strained super-lattices, superconducting films, 
and multi-layered complex structures^. Alternatively, pulsed deposition can be realized by 
chopping the flux of a continuous source with a rotating shutter 3 } independent of the exper- 
imental process we call all these time-dependent depositions PLD. Compared to ordinary 
MBE, those surfaces grown by PLD may exhibit a different surface morphology 4 . This ob- 
servation triggered several theoretical studies. Focusing on low energies there were numerical 
results that PLD crosses over to MBE 5 " 8 and that the nucleation density is characterized 
by unusual logarithmic scaling laws^, which can be motivated in terms of local scaling in- 
variance with continuously varying exponents^. Moreover, it was shown that the influence 
of a strong Ehrlich-Schwoebl barrier in PLD may lead to a smoother surface compared to 
MBE 9 . A mechanism for film in PLD is described in and review for more realistic models 
which connect with experiment | 12 | 13 [ The mean-field approach is a convenient approxima- 
tion but its limitation is described in 1 14 , 15 1 in d = 1 and in 1 13 1 16 1 for d = 2. So far theoretical 
studies of pulsed deposition have been based mostly on numerical simulations. The aim of 
the present work is to suggest a model for pulsed deposition which is expected to be valid 
in the submonolayer regime, before coalescence, and in d — 1,2 (i.e the aggregation of two 
or more islands and an adatom is neglected, only binary reactions are considered). To this 
aim we use a mean field approach for PLD found in ^ but we use the analytical rate of the 
generalized mean field in d = 1, 2P^ and we found the exact asymptotic solution of the time 
evolution of the total density of immobile islands. 

Although the modification in the rate equation is rather small, we demonstrate that it 
changes the entire scaling behavior already on the mean field level. The asymptotic solution 
for the immobile island density is t 0,5 for d = 1,2 whereas for MBE the solution is t 1 ^ 3 for 
d = 2 and t ' 25 for d = 1^. This confirms the validity of the generalized Smoluchowski 
rate for d = 1 not only for systems without a source or with a constant source but also for 
time-dependent input of particles^. 
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II. EQUATIONS FOR PLD 



We consider the following aggregation model. We assume: 

1. Brownian motion for monomers (with diffusion constant D). The islands of mass k > 1 
are immobile (i.e. Dk = for k > 1). 

2. Irreversible aggregation, i.e. when an adatom N\ contacts and thus aggregates irre- 
versibly and eventually forms an island with a larger mass +1. 

In the first approximation we assume that the island are point-like (e.g the effect of their 
lateral dimension is negligible in comparison to the effect of diffusion). This is a reasonable 
approximation up to coverage Ft such that the ratio of the average island-size to the average 
island-distance is much less than one. Hence we assume that the effective radius is r« = 1 
with a convenient choice of unit of length. Other parameters such as the capture number 
of incident monomers k s can be also introduced The model like all the Smoluchowski 
models, are valid for small time. The introduction of other parameters (r^ = i 1//d/ ,fc s ..etc) 
could improve the discrepancy we will find in the Monte-Carlo simulations in d = 2. Of 
course for more larger time scale such parameters would be relevant. But for small time our 
model, with such approximation exhibits a scaling regime for PLD. 

In our model we disregard the spatial density n^ir, t) of an island of mass k for a given posi- 
tion r and time t and compute the spatial average island density iV& = -jj^ J d^rn(v, t) =< 
n(r,t) > r . The time evolution of the average adatom density N± and the immobile island 
density N = ^2 k ^ are obtained by the generalized Smoluchowski approach. The original 
model of Smoluchowski is valid for d > 3 and, with a logarithmic correction, for d = 2. It 
is possible to use the known Smoluchowski model for d — 1, but the solutions can exhibit 
big discrepancies. The approach found in ^ generalizes the Smoluchowski model with rates 
,and the approach in ^ generalizes for the size distributions. These rates can be separated 
into two rates, the first is the the mean field rate (found in the Smoluchowski model), and 
the second is called the correlation rate. The time evolution in d = 2 is (we consider the 
reaction rate relevant for a small time scale, i.e. the mean-field rate ^ the similar equations 
were found inP^): 





1 — — 



DN 1 (2N 1 + N) +F(t) 



(1) 



N = DNl, 
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and for a very anisotropic surface, when the diffusion of adatoms is in one direction, we 
consider the d — 1 case, hence the rate is the correlation rate l 18 l 19 l(the correlation rate is the 
rate valid for more larger time than the mean-field rate, i.e. in d — 1): 

jVi = -D(N + JVi)iVi (2iVi + JV) + F(t) (2) 
JV = D(N + JVi)JVf . 

We obtain the correlation rate by multiplying the mean-field rate with the average radius 
M Q = N + Ni. The physical interpretation of the average radius is the effective surface in 
interaction with an arbitrary particle. An identical approach for MBE in d = 1 is described 
irJHl which are in contradiction with at least for the small time regime JVi <C JV. 

We assume a time dependent deposition of particles F(t) . In EDESl they consider a 
chopped deposition with two time scale, the time span of the deposition of adatoms and the 
time between the pulses. Considering that, contrary to the model proposed in the time 
between two pulses is much larger than the time between the depositions of adatoms during 
the pulses. In the limit of very short pulses the flux is of the form: 

oo 

F(t) = lY / 6(t-t k ), 

k=0 

where / is the pulse intensity defined as the density of adatoms deposited per pulse. The 
index k enumerates the pulses. For simplicity we assume that the pulses are separated by 
a constant time interval T, i.e. 

t k = kT. 

In order to find the density of adatoms N% and immobile islands iV we approximate the 
equation first for small time t and then large time t. For small time t there are few islands, 
i.e Ni ^> N, the equation reads 

oo 

JVi = 8(t - Tk) - 2DNt d (3) 

fc=0 

JV = DNt d - 

After a few pulses there are more immobile islands than adatoms. We perform the 
approximation JVi JV and the equation is: 

oo 

JVi = I 6{t - Tk) - DNiN 3 - 4 (4) 

k=0 

JV = DNlN 2 - d . 
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In the following, we will solve theses pairs of equations ([3]) and (|4]) and compare them 
with Monte-Carlo simulations. We notice that in MBE we have two parameters F, D and 
for PLD three parameters T,I,D and two unknown densities N,N%. For MBE, after a 
rescaling we have a scale free equation and for PLD the equations depend on a parameter 
which distinguishes MBE from PLD. The rate equation for d = 2 with the correlation rate 
is approximately^ the rate equation for d = 1 : 

iVi = -D(N + JVi)iVi (2iVi + N) + F(t) 
N = D(N + N x )Nf, 

we will not consider this case separately because the solutions are the same as the case d — 1. 



III. SOLUTION OF THE RATE EQUATION 

The asymptotic solution for PLD is less simple than for MBE. In order to find the time 
evolution, let us consider the temporal evolution of the adatom density between two pulses 
for t G (kT,(k + 1)T). N x is a quickly varying function in comparison to N. We then 
assume that N is constant between two pulses (i.e. for kT < t < (k + 1)T) and we perform 
an adiabatic approximation^. We solve the large time equation Q. For an arbitrary N 
we have the formal solution between the pulses(H.Hinrichsen proposed another method for 
solving these equations and found the exponent 0.5 for PLD): 

Nt = A k exp(- J D f N 3 - d {u)du) = A k exp{-D{t - kT)N 3 - d (t)), 

JkT 

where A k is the amplitude between the pulses. Since each pulse increases the adatom density 
by /, we have \im e ^ (Ni((k + 1)T + e) — Ni((k + 1)T — e)) = J, hence the amplitude A k 
follows the recurrence relation: 

A k+l = A k eM-DN 3 kT d T) + I, (5) 

where N k ? = N(kT). The amplitude is: 

1 _ e -DN 3 k ~ d kT 

A k = I 

1 —DN T 

The adatom density reads 

l-expi-DN^T) ,K k " 



Let us now turn to the temporal evolution of the density N(t) on time scales extending over 
many pulses. We consider a time scale which is large in comparison with the time T between 
the pulses, hence N is almost constant between the pulse. In order to compute the solution 
on such a large scale we will use the adiabatic approximation: 



N(t) 



N(k+i)T — N kT 
f 

f£ +1)T N(u)du 
T 

f£ +1)T DN{u) 2 - d N x {ufdu 
T 

1 - exp(-DN%- d Tk) \ 2 (1 - exp(-2 J DA^ fc 3 - d T)) 



1 - exp(-DN*- d T) J 2TN 
Two asymptotic behaviors can be observed: For a fixed N, the large time MBE, when 
Tk — t — > +oo with TN^ d — > we have N = DN lt dT2 , and the large time PLD, when 
T — > +oo we have N = 



Then the asymptotic solution for the total immobile island density, for large time, reads 

1. small T and large t the MBE regime for large time N = ^ ( 5 ~^) F t y/(5-d) 

2. large T and large t the PLD regime N = (^) 1 ^ 2 . 

For a given T the time t c defined by the equation DN(t c ) 3 ~ d T ~ 1, using one of the last 
equations for N then 

t c ~ (6) 

distinguishes the MBE regime from the PLD regime, i.e. for t c ^> T^- d the PLD regime 
dominates. We notice that, although all the equations depend on the dimension d, the PLD 
asymptotic solution is independent of the dimension, i.e. t 1//2 for d — 1,2, and the PLD 
regimes in d = 1 for an arbitrarily small time dominates (we have no MBE regime), whereas 
usually the solution is dependent on the dimensional. 
We then compute the small time regime of Eq. ([3]), hence 

Nx = -2DNf~ d for kT < t < (k + 1)T 
N = DN*~ d . 
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The formal solution for N± is(we still assume the validity of the adiabatic approximation 
e.g. the N is constant between the pulses) 

= r 5-d i i/(5-d) 
1 [ 2D(t - kT) + B k } 

where 5 fe is defined by lim e ^ (^i((A; + 1)T + e) - Ni((k + 1)T - e)) = J then: 

r 5 - d , 1/(5-4) r 5 -d -.l/(5-d) 

The asymptotic behavior for PLD is: 

Ni j- for T > 1 

Y 2D{t-kT) i -i 1^5 
L (5-4)7 i J 

and for MBE: 

Ft 

Nt — for T < 1, 

r 2D(t-fcT) i il5±d 
L (5-4)Ft ^ """J 

we have kept the first relevant correction term of the PLD and MBE regimes. Two asymp- 
totic behaviors can be observed for small time: 

1. the small time PLD N ~ D J(4 " d) t + 



\ 2D(t-kT) .1 
L (5-d)I +i J 



{ 5-d> 



2. t/ie small tame MBE when iV ~ D £ i^r i~ r- 

[ 2D(t-fcT) J ( 5=d> & fl 
L Ft(5-d) +1 J 

For the integration of the equation for the density of immobile islands we consider a time 
scale larger than T, hence 20(t _l. T) — and 2D(i _k T) — are almost constant. We recover the 

(5-d)I + 1 (5-d)Ft + 1 

small time MBE regime TV ~ f4 ^ 5 "™ 



IV. NUMERICAL SIMULATION 

In the numerical simulation of the rates following equation Eq. Q we show that the 
parameter T controls the crossover between the MBE and PLD regimes (see Fig.|TJ. We 
illustrate our approach by Kinetic Monte-Carlo simulations (KMC). In ^ a clear KMC 
realization has been presented, quite in accord with what has been presented in the present 
paper. For the KMC simulations we have a d-dimensional lattice and process as following: / 
adatoms are deposited on each pulse, between the pulses we choose randomly an adatom and 
moves it randomly in one of the 2d possible directions. The comparison between the solution 
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FIG. 1. Numerical solution of the Eq. Qwith d = 2 for PLD (T = 0.01) and MBE (with the 
same parameters F — I = D = 1, in arbitrary unit in order to illustrate the exponents). Notice 
the deviation of the PLD from MBE regime at t c = T _1 = 100 see Eq. [6] 

of the rate equation with the Monte-Carlo simulation is performed with two assumptions. 
First, we measure the nucleation density in the Monte-Carlo simulation, i.e. the number of 
island creations per-unit surface (neglecting reaction rates with more than two particles, the 
nucleation density is the same as N). Second, we consider the limit T — > oo(t kept constant), 
and in order to simulate this limit we deposit I adatoms(only on free sites on the grid) and 
when all adatoms have merged (an adatom is the neighbor of an occupied site, sticks and 
does not diffuse anymore ) to an island or an adatom, i.e. Ni(t) = then I adatoms are 
deposited (during the lapse time of deposition of / adatoms, the adatoms do not diffuse). 
In this regime the computation of the asymptotic value is simpler (see Fig. [2] and Fig.gthe 
evolution for t > 0.01 cannot be describes by our model. The author supposes that it could 
be finite-size effect or coalescence)) and we get the non-trivial scaling behavior £ 0-5 for the 
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immobile island density N. 
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FIG. 2. Kinetic Monte-carlo simulation in d = 1. Nucleation density for T = oo, e£ = 1. The best 
approximation (the dashed line given by Xmgrace) has the exponent ~ 0.518 



V. CONCLUSION AND OUTLOOKS 

In this paper we have shown that a pulsed input of particle yields a very different solution 
of the immobile cluster density. We propose a different type of crossover. We hope to 
find some experimental evidence of this scaling and crossover. The author hopes that this 
approach will lead to similar result when for instance the deposition time of particle is not 
as small in comparison to the time span between two pulses^. In this case the rescaling of 
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FIG. 3. Kinetic Monte-carlo simulation in d = 2. Nucleation density divided by t 0,5 for T = oo and 
dimension d = 2, with the density of deposited adatoms per pulse I = 10 -7 , 10~ 6 , 10~ 5 , 10~ 4 , 10 -2 
from the lowermost to the uppermost curve as function of time. 



the rates equation will lead to a set of equation with two parameters, the author expects to 
find that in this case the PLD regime dominates for large time. The PLD regime would be 
hence a regime which generalizes the MBE regime. 
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